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Abstract. Collisional excitation of the A21cm hyperfine transition is not strong enough to thermalize it in warm 
neutral ( "intercloud" ) interstellar gas, which we show by simultaneously solving the equations of ionization and 
collisional equilibrium under typical conditions. Coupling of the A21cm excitation temperature and local gas 
motions may be established by the Ly-a radiation field, but only if strong Galactic Ly-a radiation permeates 
the gas in question. The Ly-a radiation tends to impart to the gas its own characteristic temperature, which is 
determined by the range of gas motions that occur on the spatial scale of the Ly-a scattering. In general, the 
calculation of H I spin temperatures is a more difficult and interesting problem than might have been expected, 
as is any interpretation of H I spin temperature measurements. 
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1. Introduction. 

It is an article of faith among H I observers that the spin 
(excitation) temperature T sp of the ground-state A21cm 
hyperfine transition is equal to the kinetic temperature 
of the ambient gas. This belief of course arises from the 
weakness of the transition, the consequent long lifetime 
against spontaneous emission (A 2 i = 2.85 x 1CP 15 s^ 1 for 
H I), and the easily-demonstrated dominance of particle 
collisions with other H-atoms in cool gas. Much of our 
basic understanding of the phases of the ISM derives di- 
rectly from measurements of the H I spin temperature, for 
instance by comparing nearby or overlapping absorption 
and emission profiles (Dickey et al., 1978; Payne et al 
198^ ). 



But spin temperature measurements often seem to im- 
ply unphysical kinetic temperatures: not wildly so, say 
negative or infinite, but definitely in a range 1000 to 5000 
K - where the interstellar gas is unstable in multi-phase 
models ( |Wolfire et al., 1995| [McKee fc Ostriker, 1977| ) and 
therefore should be so short-lived as to be unobservable. 
One might be tempted to disparage either theory or the 
data but Davies & Cummings (1975) took a somewhat 
less doctrinaire view and showed that, for some assumed 
combinations of conditions chosen to be representative of 
gas which might produce pulsar dispersion measures, the 



gests that we remain open to the possibility that blind 
faith in the equality of the spin and kinetic temperatures 
might be misplaced. 

Calculation of the H I spin temperature turns out to be 
a remarkably complex problem, involving intimate knowl- 
edge of the ionization and phase structure of the ISM, as 
well as its topology. Here we update and expand upon 
the discussion of Davies fc Cummings (1975| ) by actually 
calculating the ionization and collisional excitation in in- 
terstellar H I regions, in the process demonstrating the 
inability of particle collisions to thermalize the A21 cm 
transit ion. But we also consider an important mechanism 
which Davies fc Cummings (1975 ) ignored, whereby the 
Ly-a radiation field threading and partly produced by the 
gas tends to impart to the A21 cm transition its own effec- 
tive temperature. Ly-a photons acquire this temperature 
while undergoing large numbers (10 7 or more) of repeated 
scatterings on many spatial scales ranging upward from 
the line-center mean free path, 0.11 (T d /100K) - 5 /iih AU 
(Td is the Doppler temperature) , and it is therefore repre- 
sentative of the motions (thermal, turbulent, etc.) which 
are established in the gas on those scales. 

That excitation by light dominates hyperfine excita- 
tion locally around individual stars is well known in the 



densitj of collision partners in warm or intercloud gas is 
simply too small to thermalize the line. Although little 
attention seems to have been paid to this warning, it sug- 



context of scattering of Solar Ly-cv radiation: Brasken & 



Kyrola, (1998), for instance, point out that Solar Ly-a ra- 
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diation dominates in local intercloud H I within 1000 AU 
of the Sun, and that the sphere of influence of an O-star 
would necessarily be much larger. In the interstellar con- 
text, we will refer to excitation of the hyperfine line by 
scattered Ly-a radiation as the Wouthuysen-Field mecha- 
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nism or WF effect following the discussions in Wouthuysen 
(19520and [Field (195$ |1959Q and we find that it opens up 
interesting possibilities for the spin temperature. At low 
density in the intercloud medium, warm neutral H I may 
(or may not) be Dopplerized by Galactic Ly-a photons but 
is certainly not thermalized by local particle collisions. In 
general, the importance which is assigned to Ly-a exci- 
tation in intercloud gas depends directly on the topology 
and large-scale structure which is assumed to apply in the 
interstellar medium. 

The plan of the current discussion is as follows. In 
Section 2 we show what is needed to calculate the exci- 
tation of the hyperfine transition. In Section 3 we lay out 
the basics of the ionization equilibrium calculations which 
are needed to determine the densities of the main colli- 
sion partners (electrons and neutral H-atoms) under given 
conditions of local gas density and gas thermal pressure. 
These calculations also serve to determine the locally- 
produced field of Ly-a photons. In Section 4 we discuss 
results of the combined ionization and particle excitation 
calculations. In Section 5 we introduce the behaviour of 
the internally- generated Ly-a radiation field, and con- 
clude that it also is not strong enough to influence the 
hyperfine excitation in neutral gas. In Section 6 we dis- 
cuss the effect of the much stronger galactic Ly-a photon 
field on the H I in various multiphase models of the ISM. 



2. Hyperfine excitation 

Discussions of the excitation of the A21cm transition gen- 



erally follow that of Field (1958) who distinguished among 
three types of excitation; radiative, by a radiation field 
around the rest wavelength of the hyperfine transition hav- 
ing a characteristic brightness temperature 7r as seen in 
the gas (Tr = 2.73K here); collisional with a total down- 
ward rate R 21 (s" 1 ) at kinetic temperature 7k ; and ra- 
diative in the Ly-a photon field, with a net downward 
rate R^, at an effective temperature 7l . Field derived 
an expression for the Ly-a excitation which was used in 
an equivalent form by Bahcall fc Ekers (1969 ) and which 
will be discussed below. 

In these terms the spin (i.e. excitation) temperature 
of the A21cm line can be expressed as 



T sp = (T R + y c + y a )/(l + y c /T K + y a /T h ) (1) 
with 

Vc = T 2l R c 21 /A 2U 
y a = T 21 R$JA 21 
and where T 2 i = hu 2 i/k = 0.068 K, A 21 = 2.85 x 

Equation 1 shows that R 2 i >> 15Tk^21 is required in 
order to thermalize the transition. 



2.1. Collisional excitation by electrons, protons, and 
neutral H-atoms 

The required rate R 21 is the sum of the downward collision 
rates over all collision partners (H-atoms, protons, and 
electrons). The downward rate constant for interactions 



with electrons is taken from Smith (1966): we note that is 
well-represented by the functional form 



log^Tk) = -9.607 + log( v /7k)exp(-(lo g r K ) 4 ' 5 /1800) 
for T K < 10 4 K and 

72i(7k>10 4 ^) = 72i(10 4 K) 

over the range considered here (the logs are base 10 
and the units of 7fi are cm 3 s _1 ) 

Smith (1966) also gives the rate coefficient for proton 
de-excitation, which is just 3.2 times larger than that for 
neutral atoms at 7k > 30 K. We used 7^ = 3.27^ but 
employed the more recent neutral atom de-excitation rates 



of Allison & Dalgarno (1969) which we fit in piece- wise 
continuous fashion, following their tabular presentation. 
Excitation by protons is unimportant because it is so much 
weaker than that by electrons at the same temperature 
and because n e > n p under all conditions. 

The results of Allison & Dalgarno (1969| ) are given up 
to 1000 K at which point they are increasing as Ik ' 33 , 
but with a negative second derivative. The neutral-atom 
rate coefficients of Smith (1966) increase somewhat faster 
for 7k > 1000 K with a slightly positive second deriva- 
tive, but are not substantially different in the mean. We 



extrapolated the rate coefficients of Allison & Dalgarno 
(1969[ ) t° higher temperature as 7k ' 33 . In this regime the 



excitation is increasingly dominated by electrons and the 
results reported here are not noticeably influenced by the 
manner of the extrapolation. 

2.2. Excitation by the Ly-a photon field 

Recoil - momentum conservation - at each scattering 
causes the field of Ly-a photons to have a slope (color 
temperature) at the line center which corresponds to the 



Doppler tempe rature of the ambient medium ( Field, 1959 



Adams, 1971a ); Ly-a photons sample the gas kinemat 
ics and then tend to impart to the hyperfine transition 
a corresponding excitation temperature. Bahcall & Ekers 
Kl969| ) and [Urbaniak fc Wolfe (198l"[ ) employed this mech- 



anism to show how the H I spin temperature could remain 
low, rendering A21cm absorption detectable, even in neu- 
tral gas near quasars. In cases where the gas kinematics 
are affected by non-thermal motion, this will be reflected 
in the spin temperature, depending on the relative impor- 
tance of collisional and Ly-a excitation. 

Field (19138) considered detailed balance in the various 
transitions between sublevels which constitute the Ly-a 
line and arrived at the following expression for y a 
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Fig. 1. Ionization equilibrium calculations in atomic gas. 
The full curves are plots of the ionization fraction n e /nn 
where nu is the number density of H-nuclei in all forms. 
Results are shown for four values of the total gas pressure 
P/k = 750, 1500, 3000, 6000 cm" 3 K and the kinetic tem- 
perature is traced across each set of curves. The shaded 
curves are plots of the electron density. Open symbols rep- 
resent the ionization fraction in warm neutral and warm 



ionized gas in the three-phase model of McKee & Ostriker 
(1977[)] The filled rectangles show the ionization fraction of 



10 



if) 



warm neutral gas in the two-phase calculations of |Wolfirc 
et al. ( [1995 ) for their standard case Nw 



cm 



y a = (cr 21 /36^)(m H /fcr D ) u - 5 A 3 (A Q /A 21 )n Q (2) 

where A a is the spontaneous emission rate (4.7 x 10 s 
s^ 1 corresponds to the mean Ly-a absorption f-valuc 
f=0.4167), n a is the number density of Ly-a photons at 
the line center and Td, the Doppler temperature, repre- 
sents the width of the Ga ussian core of the lin e. Td is 
defined in the usual way ( [Leung fc Liszt, 1976| ) by the 
expression 



2kTr>/rriH = 2kT K /m H + v tu 



rb 



leading to Td = Tk+121 K (v tMr f,/l kms 



(3) 
This 



expression accords with the definition used by Payne et al 
(19820TD = 21.86K W 2 where W is the FWHM (kms"' ) 
and differs from the earlier usage of Dickey et al. (197S| ), 
T D = 121K W 2 taking W as the HWHM The Doppler 
temperatures shown by Dickey et al. (1978 ) are actually 
too large by a factor 2*ln(2) w 1.4. 

Evaluating the constants in equation (2) we find 



Species 


a (cm d s 




b 


H n >i 


3.50 x 10" 


-12 


0.75 


H2 S 


1.54 x 10" 


-13 


0.54 a 


He 


2.36 x 10" 


-12 


0.64 


C 


4.67 x 10" 


-12 


0.60 



Rate into the 3 level only, see Section 3.2 
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5.90 x 10 n n Q T D 
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3. Ionization equilibrium in H I regions 

3.1. Calculations 

To obtain number densities of the various collision part- 
ners which determine the equilibrium of the hyperfine 
transition it is necessary to calculate the ionization equi- 
librium of a di lute gas. We duplica ted the relevant portions 
of the work of Wolfirc et al. (1995 ) which is a very compre- 
hensive study of two-phase equilibrium in the interstellar 
medium and must be regarded as the standard reference 
on this subject. That is, we consider a gas consisting of 
hydrogen, helium and carbon, ionized by soft x-rays and 
cosmic rays (and their secondary electrons) which recom- 
bines in both the gas phase and on small grains (following 
the formalism of Drainc & Sutin (1987)). This calcula- 
tion is fully described in the original reference and a de- 
tailed discussion will not be repeated here. We note that it 
uses Solar metallicity, no appreciable gas-phase depletion 
of carbon ([C]/[H]= 3.3 x 10~ 4 ), and a primary cosmic ray 
ionization rate (^ = 1.7 x 10~ 17 s _1 . The ionization frac- 
tions are overestimates by factors of 2-3 in higher-density 
gas if gas-phase depletion is important. 

Wolfire et al. (1995) parametrized their model in terms 
of a column density of hydrogen nuclei Nw which attenu- 
ates part of the soft x-ray flux; their standard model takes 
a value Nw = 10 19 cm -2 for this quantity, with factor of 
ten variations in either sense as subsidiary examples. To 
model the x-ray ionization rate, we calculated the total 
attenuation at each energy over a gas column Nw (due to 
all the important elements, not just those whose ioniza- 
tion is calculated) using the absorption cross-sections of 



Balucinska-Church & McCammon (1992) and modulated 



the free-space x-ray spectrum accordingly. The product of 
this attenuated and modulated x-ray spectrum and the 
absorption cross section is then integrated over energy to 
give the primary ionization rate for hydrogen and helium. 
We also calculated the mean photoelectron energy for use 
in determining the secondary ionization; this is typically 
40-50 eV as for cosmic-ray events. 

3.2. Recombination coefficients 

We took recombination coefficients for various species 
from the 1999 update of the UMIST reaction database. 
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Fig. 2. The spin temperature of the H I A21cm transition 
resulting from the same conditions whose ionization equi- 
librium is shown in Figure 1. Symbols represent results 
for warm gas in the two and three phase calculations of 
Wolfire et al. (1995|) and |McKee fc Ostrikcr (197% , see 
Figure 1. 

They are expressed in the form ay = ay x (Tk/SOO) - ^ 
where the coefficients are given in Table 1. Hydrogen of 
course is a special case whereby the needed recombination 
coefficient is that into the n=2 state and higher, since re- 
combinations into n=l produce ionizing photons which are 
immediately reabsorbed in the surrounding gas. A special 
case is considered below whereby we require the recombi- 
nation rate for production of Ly-a radiation. In this case, 
recombinations into the metastable 2S state must be ig- 
nored, and the effective recombination rate is diminished 
by an amount which is given in the second entry of the 
Table, taken from [Martin (198g|). 



3.3. Photoionization of carbon 

For this we take from the UMIST database the free space 
value Cc(0) = 3. Ox 10~ 10 s _1 . This can then be traced into 
a gas column as Q(Ay )= ^(0)exp(— 2A2Ay) following 
Black & Dalgarno (1977) but we will consider only the case 
Ay — 0. The major contribution to the ionization fraction 
for the warm gas considered here is generally from H and 
He, not from carbon. 



3.4. Collisional ionization of hydrogen by electrons 

We used the electron-im pact collisional ionization rates of 
Scholz fc Walters (1991 ) . The reader is cautioned to ignore 



an unfortunate 10 order of magnitude error in labelling of 
the vertical scale of their Figure 1, lest this crucial process 
be ignored. 

4. Ionization equilibrium and hyperfine excitation 
in the interstellar hydrogen 

4.1. The ionization fraction as a function of gas 
pressure and density 

In Figure 1 we show the results of our ionization calcula- 
tions for four values of the gas thermal pressure, using the 
standard value Nw = 10 19 cm" 2 . The ionization fraction 
is shown as solid lines and the electron density is shown 
shaded. The kinetic temperature is traced across the set 
of four curves of varying pressure. Clearly, the ionization 
fraction in neutral gas is controlled by the total density, 
while the electron density - which varies remarkably little 
with 7Jh until the gas ionizes fully - is influenced somewhat 
by the pressure. The soft x-ray ionization rate scales the 
overall behaviour, albeit somewhat weakly (not shown). 

The filled dark squares in Figure 1 represent the lim- 
iting cases at high and low pressure for the model of two- 
phase equilibrium from |Wolfirc et al. (1995p . Table 2 gives 
the range of pressure and kinetic temperature over which 
two phases can coexist for various values of Nw and it 
shows that two-phase equilibrium is possible over a wider 
range in Tk than is sometimes considered when interpret- 
ing spin temperature measurements. Comparison of the 
Table entries for the standard case Nw = 1 x 10 19 cm" 2 
with the locations of the filled dark symbols in the Figure 



shows that our results mimic those of Wolfire et al. (1995), 
as intended. Both sets of calculations produce the same 
ionization at any given density and temperature. 

In Figure 1, two unfilled s ymbols represent the results 
of McKee fc Ostrikcr (1977 ) for the warm neutral and 
warm ionized gas (both at 8000K). The pressure in the 
McKee-Ostriker model is P/k = 3700 K cm" 3 but the ion- 
ization fraction in three-phase warm neutral and ionized 
gas is considerably highe r than that in our calculations. 
McKcc fc Ostrikcr (1977| ) explain in the caption to their 
Figure 1 that the warm neutral medium can only be pro- 
duced in appreciable quantity at the assumed pressure by 
assuming a higher than average value of the soft x-ray flux. 
In the context of our modelling, the soft x-ray flux would 
have to be increased by a factor of about 50 in order to 
give the same total ionization which occurs in warm, neu- 
tral gas in the McKee-Ostriker model. It is unfortunate 
that the basic parameters of the three-phase model have 
not been revised in so long. 

4.2. Evaluation ofT sp 

Figure 2 shows that for kinetic temperatures above 2000 
K, and especially at 8000 K, the particle excitation rate 
is too small to thermalize the A 21cm line in any model 
of the ISM, including the McKee-Ostriker three-phase pic- 
ture. The extreme conditions assumed for the three- phase 
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Table 2. Two-phase spin and kinetic temperatures 



Model 


P/k (cm d K) 


T K (K) 


T SP (K) 


Nw=10 18 cm" 2 


1600 - 9990 


9200 - 5600 


3525 - 4620 


Nw=10 19 cm" 2 


990 - 3600 


8700 - 5500 


1800 - 3260 


Nw=10 20 cm" 2 


610 - 1500 


8200 - 4900 


1035 - 2020 


Three-Phase 


3700 


8000 


5400 



model result in T sp « 5400 K in warm neutral gas with 
Tk= 8000 K. In the two-phase model, the spin temper- 
ature of warm gas would be smaller, 1800 - 3200 K for 
the standard conditions where 7k= 8700 - 5500 K (the 
models at higher Xk have lower pressure and lower T sp ). 
Table 2 shows the range of spin temperature which results 
from various conditions of two- and three-phase equilib- 
rium. Figure 3 displays the y-factor for collisional excita- 
tion corresponding to the conditions shown in Figures 1 
and 2. In the warm gas, y c /T^ « 1, while j/c/Tk >> 1 is 
the relevant condition for thermalization. 

Physically, the subthermal excitation comes about be- 
cause the condition for thermalization becomes ever more 
stringent at higher kinetic temperature, while the neutral 
particle excitation rate coefficient only increases as Tk 



(Allison & Dalgarno, 1969). The neutral particle density 
declines as 1 /Tk at constant pressure until the gas begins 
to ionize, at which point new particles enter the gas and 
the decline is even more rapid. Once the gas ionizes even 
slightly (a few percent), excitation by electrons begins to 
dominate and the remaining neutral atoms will be more 
strongly excited. But the transition from neutral to ion- 
ized gas (which increasingly precludes observation of the 
A21cm line) occurs at too low a density to preserve the 
thermalization of the line even when a very large soft x- 
ray flux is asserted, as in the model of McKee fc Ostrik er 
(1977^B 

We conclude that in the context of these calculations, 
one simply cannot expect to observe spin temperatures 
which are directly indicative of the 8000 K kinetic tem- 
perature typically assumed (by observers) for warm gas. 
Instead, seemingly unphysical, intermediate values will ap- 
pear, even when the contribution of warm gas can be iso- 
lated. This situation comes about mostly because the ex- 
citation is sub-thermal, but also because multi-phase equi- 
lbrium is not entirely oblivious of local conditions; the ki- 
netic temperature of warm neutral gas can vary between 
5000 and 9000 K in two-phase models depending on the 
pressure, elemental abundances, etc. Moreover, Wolfire 



et al. (1995) note that the time to reach thermal equilib- 
rium in warm gas is too long to ensure that a steady-state 
actually occurs; rather, two-phase equilibrium is an ideal, 
toward which the ISM may tend. If seemingly forbidden 
spin temperatures 1000-5000 K are actually measured in 
the A21cm line, it is unfair and improper to impugn the 



1 Perhaps this model might also be referred to as the Lake 
Woebegon solution, since it seems to assume that the soft xray 
field is everywhere brighter than average. 



multi-phase models on this account, because they predict 
them. 

Finally, we note that the need to consider a wide range 
of spin temperatures in warm gas makes it harder to dis- 
cern T sp in cool gas when both are seen blended along 
the line of sight. Of course, lines of sight in which the 
contribution of cool gas is not strongly blended are al- 
most unknown. The presence of lower T sp in the intercloud 
medium would result in lower T sp in line blends, perhaps 
producing another range of seemingly unphysical T sp be- 
tween w 1000 K and the range e xpected for cool gas 
alone, 40-200 K (|Wolfirc ct al., 1995|) . 



5. Excitation by the Ly-a photon field 

In the previous section we showed that the A21cm hyper- 
fine line will not be thermalized by particle collisions in 
warm weakly ionized gas. But there is another means of 
excitation available, the Wouthuysen-Field effect, which 
we now discuss. 



5.1. Ly-a scattering 

Our basic understanding of the behaviour of Ly-a photons 
in very optically thick media dates from the discussion of 



Field (1959) on relaxation in Doppler broadening and the 



numerical work of Adams (1972) on scattering including 
strong damping. The latter showed that the mean number 
of scatterings (N) of a photon before escape from a loss- 
less medium never exhibits the behaviour expected of a 
pure one-dimensional random walk. Taking to as the line- 
center optical depth at the midplane of a semi-infinite slab 
viewed along the outward normal, Adams fo und instead 
that, at suffic iently high optical depth (see Hummer fc 
Kunasz (1980|) ), (N) « tq. This seems to have been un- 



expected, even though identical behaviour is discussed by 
Field (1959 ), citing ^anstra (1949 ), following essentially 
the same line of reasoning. 

Adams (1972) explained this result heuristically by 
noting another aspect of the simulations, namely that 
photons escape the medium only in the far line wings 
of the Voigt profile </>(a:), where x is a normalized fre- 
quency displacement from line center just below. They 
escape at \x*\ = |(^» — v )\/Ai>d w (aro) 1 / 3 where 
Ais D = K/c) y/2ln(2)kT D Jm H is the Doppler HWHM 
and a is the damping constant measured in units of 
Ave>. Although most scattering occurs nearer the Doppler 
core, photons move appreciably in space only when 
they scatter at frequencies where the medium is less 
opaque. The photons escape at a frequency where the 
optical depth is still well above unity; |x*|/|av = i| = 
O^iTo/W^K^^iNH/W 20 cm" 2 )" 1 / 6 where ±x T=1 are 
the frequency shifts at which the optical depth drops to 
unity for a pure damping profile. 

Given that the rms Dopper shift in any scattering is 
|A^£>| (the mean shift is in the core and — l/\x\ in the 
far wings), Adams reasoned that, occasionally, a photon 
would experience one series of scatterings into the far 
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Fig. 3. The y-parameter for collisional (solid curves) and 
Ly-a excitation (shaded and chained curves) for a medium 
of assumed size 2R = 3 pc. The same 4 thermal pressures 
are employed here as in Figures 1 and 2. The shaded curves 
marked by "R=1.5pc" represent the correction factors, ac- 
counting for radiative transfer in lossy media, which have 
been applied to the Ly-a excitation following the discus- 
sion in Section 5.2; see also the following Figure. The 
horizontal line at lower right shows where Nh = 2iihR 
cm" 



wings carrying it at once frequency units from the 
line center and entirely out of the scattering medium. If 
a photon must travel a distance tq to escape, x* should 
satisfy the condition x*/<j)(x*) = tq since the scattering 
length is proportional to l/4>(x). This recovers the be- 
haviour seen in the simulations since 4>{x] oc a/x 2 in the 
damping wings. 

Shortly after this work, the tra nsfer equation was ef- 
fectively solved in closed form by Harrington (1973 ) in 
the case of very large optical depths (ar ) > 1 The 
behaviour seen in the simulations of Adams (1972 ) was 
recovered analytically in the change of variables needed to 
effect the solution: to something which varies as x 3 in the 



line wings. The numerical work was extended by Adams 
(1975[)]and Hummer fc Kunasz (198C ), and the analytic 



work by |Neufeld (1990| , |l99l| ). There, it was shown that 
the mean distance travelled by a Ly-a photon is actually 
(ar ) 1/3 r . Because a = 0.0425/T D °' 5 is typically a small 
number, the cube root factor is of order 10-100. This de- 
scription is correct to factors of order unity, as discussed 
by Hummer & Kunasz (1980). 

For reference, we note that the line center optical 
depth T = 5.90 x 10~ 12 7V H /Tb°- 5 is proportional to 



A a /Ai>]j, which is the damping constant a except for 
a numerical factor. Thus gtq = 1.39 x 10~ 10 a 2 iVH = 
2.51 x 10- 13 7Vh/Td- If T d ~ T K , the latter ratio of column 
density to temperature is familiar from the expression for 
the integrated optical depth of the hyperfine line. Thus 
we come full circle and find (but do not intend to claim 
a causal connection!) that the crucial parameter for Ly-a 
excitation of the A21cm line is proportional to the A21cm 
optical depth. 

5.2. Internal sources of Ly-a photons 

We included three locally generated sources of Ly-a pho- 
tons in H I gas: electron impact excitation of the 2 P level 



using the rate constant calculated by 3cholz & Walters 
|(199lD ; electron- proton recombinations into levels n > 2 
(except 2 S) for which the relevant coefficients are summa- 
rized in Table 1 and for which our rates are very slight 
upper limits; and excitation by the same cosmic- and x- 
ray fluxes which ionize the gas. For the latter, we follow 



Draine & Salpeter (1978) and assume that 1.5 Ly-a pho- 
tons are generated per primary ionization event. 

In general, electron-impact excitation is unimportant 
except above 10 K and cosmic-ray and x-ray generated 
photons dominate in the denser gas. 

5.3. The locally-generated Ly-a photon field in a dusty 
H I cloud 

The volume rate at which Ly-a photons are created, n a , 
is straightforward to specify. But for the excitation rate 
we require in Field's formulation the volume density of 
photons n Q = n a x t, where t is a characteristic retention 
time (the expression could just as well have been c ast in 
terms of the local mean intensity). In fact, following Field 
1958| ), the relevant time is R/c where 2R is the character- 



istic size of the medium. Although this seems naive, it is 
supported by the analytic solutions of Harrington (1973) 
which we now discuss: the actual distance travelled by a 
photon, (aroY^TQ, is associated with a broadening of the 
profile by an amount (aTo) 1 / 3 , leaving the line-center Ly-a 
intensity unchanged. 

Harrington (1973) presents an analytic solution for the 
case of a uniform slab with a uniform internal source 
of Ly-a radiation h a such that 1 photon is created per 
unit area per unit time looking down through a slab of 
half-thickness R. To first approximation, the mean inten- 
sity at the slab midplane, assuming a lossless medium, 
is J(t = 0,x) = (\/6/27r 3 )exp(-1.34|x| 3 /(aT )). To con- 
vert this normalized solution to "real" units, the density 
of photons at the line center is just Air/c x 2R x J(0, 0) = 
(4\/6/7r 2 ) x (R/c) = 0.993i?/c= / x R/c. 

A slightly more accurate expression for the photon field 
is given by Harrington (1973 ) as a series which includes the 
possibility of accounting for losses by continuous absorp- 



2 The discussion of the previous two paragraphs neglects in- 
ternal absorption, which is discussed in detail in Section 5.3 
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Fig. 4. A scaling factor for the mean intensity inside 
a lossy slab harboring a uniform source of Ly-a radia- 
tion, following equations 27 and 41 of Harrington (1973] ); 
see Section 5.3 here. Values are shown at the midplane 
(t/to = 0) and two-thirds of the way from the midplane 
to the slab surface. Relevant values for warm and cool neu- 
tral gas are indicated. These calculations are the origin of 
the correction factor displayed in the prior Figure. 



tion in the medium. To calculate the loss it is necessary 
to evaluate the term (Vo/2)ero where e is the (assumed 
small) probability of absorption per scattering (recall that 
there are on average To scatterings). The mean distance 
travelled per scattering, measured in units of the distance 
necessary to traverse one optical depth at the line center, 
is (ar ) 1/3 flAdams, 1975] ; |Neufeld, 1990| ); the physical dis- 
tance corresponding to one optical depth at the line center 
is found by equating Ar = 5.90 x 10 _12 n H x AI/^/Td = 1. 



From this it follows that etq — {aTo) 1 / 3 N-na gr abs where 



lbs 



• gr 



= 1.27 X 10" 



gr 

cm -2 has been calculated re- 



cently by Weingartner fc Draine (2001 ) (their value for 
the albedo is 0.32). Several c ases of lossy media are treated 
at length by Ncufcld (1990] ) (analytically) and by Adams 
(1971a|) (numerically). The formalism discussed here re- 
produces the numerical results in the latter reference quite 
well but the grain model of Weingartner fc Draine (2001 ) 
is about 30% more absorptive than the worst case consid- 
ered by Adams. 

In Figure 4 we show the scale factor f multiplying R/c 
as a function of ero. For a lossless medium, the exact an- 
alytic solution yields f=0.909 at the line center and slab 
midplane(vs 0.993). Losses are relatively unimportant in 
the intercloud gas over columns which are only a few pc 
long, i.e. those which might be associated with individual 
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Fig. 5. As in Figure 2 when a specified fraction F Q of 
the galactic flux from early-type stars threads the gas and 
T\ = Tk- The panel with F Q = repeats behaviour already 
shown in Figure 2. 



cloud envelopes (the case where the size of the medium 
is comparable to the galactic scale height is considered in 
the next Section). But the attenuation of Ly-a radiation 
will be significant in any cool gas which absorbs strongly 
at A21cm. 

5.4. Efficacy of the locally-generated Ly-a photon field 

Figure 3 is a behind-the-scenes look at the excitation of 
the hyperfine line as it occurs under the conditions used 
to produce Figures 1-2; it shows the y- values for colli- 
sions and Ly-a excitation. For the purposes of the calcu- 
lation, we took a characteristic size R = 1.5pc. Ly-a ex- 
citation would compete strongly with particle excitation 
in cooler gas were it not for the effects of dust absorption 
and the decrease in ionization fraction due to recombina- 
tion on grains. In warmer gas, the small physical size of 
the medium assumed for the calculation limits the effec- 
tiveness of Ly-a excitation except where the gas ionizes 
fully. 

We conclude that there is no locally-generated source 
of excitation within individual clouds which will suffice to 
thcrmalize the spin temperature of warm H I, although 
the small attenuation in warm gas means that Ly-a exci- 
tation can become important if galaxian scales are consid- 
ered. Because the Galaxy is a copious producer of Ly-a 
photons, the galactic Ly-a radiation field will dominate 
excitation of the A21cm line under some assumed condi- 
tions, as discussed below. 
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Fig. 6. Variation of the Ly-a radiation field inside a slab 
of warm H I harboring a unit plane source 5(t) at the 
midplane, plotted against distance from the slab center 
t/to; see Section 6.1 of the text. Results for a lossless slab 
are shown, and for Nh = 1, 2, and 3 x 10 20 cm -2 . 



6. Ly-a excitation by Galactic radiation 

The space-averaged volume creation rate of Ly-a photons 
near the galactic plane dwarfs that generated in situ in H I 
gas by large factors. The galactic density of Ly-a photons 
is not known with certainty but estimates are enough to 
force consideration of its importance in exciting the neu- 
tral H I in warm media. 

Vacca et al. (1996) have shown that the areal pro- 
duction rate of ionizing photons from 429 O- and early- 
B stars within 2.5 kpc of the Sun is * = 3.7 x 10 7 
photons cnr 2 s _1 . Each ionizing photon should pro- 
duce of order one Ly-a photon. Thus, over a cylinder of 
height h = ±50 pc = 100 pc about the Galactic mid- 
plane, the mean volume production rate is n a w ^o/h = 
1.2 x 10~ 13 cm~ 3 s _1 . This is far higher than that gen- 
erated solely by recombination in the warm neutral gas 
we modelled: the combination of a long path and higher 
volume production rate would result in a vastly increased 
line excitation by Ly-a photons. ^ 



This value of typ is three times higher than that of 



Drainc 



&i Salpeter (1978) who added the separate contributions of H 



II regions, recombination in diffuse hydrogen, and production 
in supernova rema nts (any such e xtra contributions should in 
fact be included). Adams (1971b |) cites an observational limit 



on the galactic Ly-a flux due to [Thomas fc Krassa (1971 ) but 
this limit is consistent with large values of given recent 
estimates of the continuous absorption by dust. 



In Figure 5 we show the behaviour of T sp when a spec- 
ified fraction F a of the galactic Ly-a flux threads the gas 
assuming Nw= 10 19 cm -2 and 7l= 7k • Behaviour in the 
upper left panel of Figure 5 is identical to that in Figure 
2. Figure 5 shows that even a very small fraction of the 
galactic Ly-a radiation suffices to dominate the excitation. 

Apparently, there is an abundant source of non- 
collisional excitation, if only it can be used. 

6.1. A two-phase ISM 

We consider first an idealized two-phase interstellar 
medium which fills all space, and in which nearly all of 
the volume is occupied by warm i ntercloud gas. F ollowing 
the discussion of this problem by Adams (1971a) we ide- 



alize the thin layer of young stars to be a plane source 
embedded in the much thicker, warm H I layer, and apply 
an analytic radiative transfer solution using the new value 
of the grain a bsorption cross-section from Weingartner & 
Drainc (2001 ) as mentioned earlier. At first, we neglect 
the fact that most photons which enter cold H I clouds 
will be absorbed by dust. 

Figure 6 shows our evaluation of the analytic solution 
( Harrington, 197S| ; Ncufcld, 1990) to this idealized prob- 
lem, performed for a gas with 7k= 7b= 7200 K and var- 
ious values of the column density in the H I layer (which 
determines the loss due to continuous absorption by dust). 
Although the total column density of H I looking out from 



the midplane is about 3 x 10 cm 2 flLiszt, 1983 ), that of 



the intercloud gas is somewhat less, of order half. Figure 
6 shows that over most of the volume of this idealized 
model, the Ly-a radiation field threading the intercloud 
gas is certainly not less than 1% of ^o- In this case, the 
spin temperature in the intercloud gas at all but the high- 
est galactic z-heights would be close to equilibrium with 
the temperature of the ambient galactic Ly-a radiation - 
7k, 7b or what have you, depending on how the photon 
field relaxes in its interactions with the gas. 

Cool neutral gas packets dispersed throughout such a 
scattering medium would clearly be threaded by a sub- 
stantial fraction of the galactic flux as well, raising the in- 
teresting possibility that the spin temperatures measured 
in strongly absorbing gas are also influenced by Ly-a ra- 
diation. A calculation analogous to that used to produce 
Figure 5 shows that Ly-a excitation dominates if more 
than 0.1% of the galactic radiation field threads cool gas. 
In this case, T sp = 7b rather than T sp = 7k and typical 
turbulent velocity contributions to Td will put a 30 - 120 
K floor on measurements of T sp (see Section 6.3) if the 
turbulence is established on the scale of the Ly-a scatter- 
ing. 



6.2. A multi-phase ISM 



Neufeld (1991) pointed out that a multiphase ISM struc- 
ture may foster escape of Ly-a radiation and that indi- 
vidual neutral gas clouds - warm or cool - may not be 
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Fig. 7. Line-center Ly-a mean intensity inside a lossless 
slab in the presence of a half-unit monochromatic flux 
6(v — Vj) /2 incident on e ach face following the analytic so- 
lution of Neufeld (199C ) . Two combinations of slab thick- 
ness and Doppler temperature are shown, representing 
cool and warm neutral gas. For each of these, the mean 
intensity is shown at the center and 90% of the way from 
the mid-plane to the slab surface. Over most of the in- 
terior, the behaviour more closely resembles that at the 
midplane. 



threaded by much of the ambient Ly-a flux. If most of the 
volume of the ISM is in a contiguous hot phase having 
negligible neutral hydrogen or Ly-a absorption, scattered 
radiation escapes relatively easily even when neutral gas 
clouds have a very high surface covering factor. Photons 
scatter off individual clouds as if they were particles of 
a gas having a Doppler temperature corresponding to the 
cloud-cloud velocity dispersion, 6 kms" 1 , appreciably nar- 
rowing the interstellar line profile if the photons are scat- 
tered sufficiently often before escape. The crucial issue is 
the fraction and connectedness of the ISM occupied by 
very hot gas in which hydrogen is totally ionized, as op- 
posed even to 6000-10000K H I intercloud gas which will 
still present a sizable opacity to Ly-a radiation. 

The fraction of incident photons which will be reflected 
back from the outer boundary of an isolated, discrete cloud 
is of order 1 — 1/(to0(x)) where x is the normalized fre- 
q ucncy shift o f the incident radiation from line center 
see Section 5.1 here. An analytic solu- 



flNcufeld, 1990|) : 
tion to the prob lem of radiatio n incident on a lossless slab 
is presented by Neufeld (1990| ) in his equations 2.19-2.22; 
calculation is considerably simplified when Neufeld's equa- 



tion 2.22 1^| is recast as 

F(w, y) = g x ln(l — 2 exp(— ny) cos(ttw) + cxp(— 2ny)), 

g = V6/(16tt 2 ). 

In Figure 7 we show the mean intensity at the line 
center inside lossless slabs of warm and cool gas as the 
frequency shift of the attempted injection (expressed in 
kms -1 ) varies; neglect of loss is appropriate for warm gas 
on the scale of individual clouds. Any cool gas in the ISM 
is either embedded in a strongly-scattering intercloud gas 
or enrobed in warmer H I and its internal radiation field 
could be much larger than that shown here; the reflectivity 
calculations for bare cool gas are included only to illustrate 
parameter sensitivities. 

If the ambient Ly-a photon field is as narrow as the 
6 kms -1 cloud-cloud velocity dispersion after multiple re- 
flections, little will penetrate cool gas and essentially none 
will enter the warm H I. But if the Ly-a photon field is 
appreciable at relative velocities of 30-100 kms -1 with re- 
spect to the line center, the internal radiation field in warm 
neutral gas could be as high as 0.1-1% of that which is in- 
cident. The interstellar Ly-a field would be quite broad 
in any model where the galactic layer harbors apprecia- 
ble amounts of distributed warm H I (Adams, 1971b) and 
Ly-a photons will be injected into the ISM in the local 
line wings when they escape an H II region. But the am- 
bient Ly-a radiation field and its Doppler temperature 
are not easily specified in extreme three-phase models in 
which no Ly-q scattering occurs between individual clouds 
flNeufeld, 1991]) . 



6.3. Non-thermal Doppler excitation 

For cool gas, where Ly-a excitation would be important if 
galactic radiation enters via the intercloud or warm neu- 
tral gas which enrobes it, HI absorption features typically 
show a noticeable difference between T sp and Td, imply- 
ing that excitation of the hyperfine transition is not dom- 
inated by Ly-a excitation at a temperature other than 
T K - 

In Figure 8, we show the turbulent velocity contribu- 



tions derived from the measurements of Payne et al. (1982 



1983 ) where, for the stronger features (r > 0.04) and ex- 
cepting a few cases with very large quoted errors in T sp , 
we have derived values of Vturb according to equation 3. 
In those few cases where Td < T sp in the data, the value 
is shown negative. The error estimates are ours, derived 
from analysis of the errors quoted in the original data. 

Clearly, most of the data are consistent with an 
added microturbulcnt contribution of order Vturb = 0.5- 
1.0 kms -1 , equivalent to adding a 30 - 120 K turbulent 
contribution in quadrature with 7k to form Td- Of course, 
if Ly-a excitation is actually operating at a Doppler tem- 
perature above 7k, the measured spin temperatures are 



4 Note that the quantity plotted in Figure 2 of |Neufeld (199d| ) 

is (aro) 1//3 J(r , x), not (aro) J{to, x). 
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Fig. 8. Turbulent velocity contributions derived from the 
spin and Doppler-temperature measurements of Payne 
et al. fll982| , |1983| ) using (T D = T sp +121 v t . 
tures having tq > 0.04. 



urb 2 ) for fea- 



already somewhat inflated, and the turbulent velocities 
extracted from the comparison between T sp and Td are 
only lower limits. Statistically, (T sp /Tb) = 0.56 ± 0.26 in 
this dataset. 

Recent studies of stray-radiation-corrected Galactic 
emission profiles find evidence for three kinematic compo- 
nents having lincwidths of 5, 13, and 31 kms -1 (Verschuur 
& Perftt, 1999). These correspond respectively to individ- 
ual clouds (or to Tu= 550 K; see Takakubo (1967) and 
references therein), to the 6 kms -1 cloud-cloud velocity 
dispersion (or Td= 3700 K), and to an anonymous com- 
ponent with Td = 21000 K or a combination of Tk = 
8000K, vturb = 10.4 kms -1 . Suggestions of a pure inter- 
cloud component corresponding to Tk = 8000K, v tU rb = 
kms -1 have a somewhat checkered history, but seem now 
not to be clearly evident in the data. This is perhaps not 
too surprising given the panoply of phenomena which may 
influence an observed linewidth on galactic scales. 

7. Summary 

The major conclusions of this work may be summarized 
as follows: 

Particle excitation does not thermalize warm (inter- 
cloud) H I in cither the two- or three-phase models of the 
ISM. In the absence of other sources of excitation, one 
would expect to see T sp sw 1000 — 5000 K, from regions 
where Tk ~ 5000 — 10000 K (note that regions of higher 

typically are at lower pressure in two-phase equilib- 



rium and produce lower T sp ), even when the contribution 
of warm gas can be isolated. The presence of a wide range 
of T sp in warm gas complicates the derivation ofT sp in cool 
gas, which is nearly always seen blended with warm gas 
along the line of sight. 

In the two-phase model, at all but the highest galactic 
z-heights, hyperfine excitation in warm, neutral intercloud 
gas will generally be dominated by the scattering of galac- 
tic Ly-a photons produced by the ensemble of OB stars 
and supernova remnants, in which case the spin temper- 
ature will correspond to whatever Doppler temperature 
the Ly-a photon field acquires while interacting with the 
interstellar gas. Alternatively - recalling comments in the 
Introduction - excitation may be dominated locally by the 
Ly-a radiation originating in or near particular objects. 

In the three-phase model the ability of the galactic 
Ly-a radiation field to penetrate individual H I clouds is 
somewhat problematic. 

In any multi-phase model where the warm neutral gas 
is threaded by a substantial galactic Ly-a flux, excitation 
in the cool H I will be influenced by the Ly-a photon field 
also. 

The excitation in cool H I clouds seems not to be to- 
tally dominated by Ly-a excitation at a temperature other 
than Tk, judging from the ubiquity of an apparent turbu- 
lent velocity contribution to the H I absorption linewidth, 
i.e. (T sp /T D ) « 0.56. 

The scattering calculations which will tell us the color 
temperature of the interstellar Ly-a radiation field remain 
to be performed and may require more knowledge of the 
structure of the ISM than is currently available to us. 
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